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We study the dynamics of Dirac-Born-Infeld (DBI) dark energy interacting with dark matter. 
The DBI dark energy model considered here has a scalar field with a non-standard kinetic energy 
term, and has potential and brane tension that are power-law functions. The new feature considered 
here is an interaction between the DBI dark energy and dark matter through a phenomenological 
interaction between the DBI scalar field and the dark matter fluid. We analyze two different types 
of interactions between the DBI scalar field and the dark matter fluid. In particular we study the 
phase space diagrams of and look for critical points of the phase space that are both stable and lead 
to accelerated, late-time expansion. In general we find that the interaction between the two dark 
components does not appear to give rise to late time accelerated expansion. However, the interaction 
can make the critical points in the phase space of the system stable. Whether such stabilization 
occurs or not depends on the form of the interaction between the two dark components. 

PACS numbers: 04.60. Cf; 98.80.Es 
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I. INTRODUCTION 



In recent years much of the effort in theoretical physics 
has gone into the study of the observed present acceler- 
ated expansion of universe first reported in P, |3| through 
observational data from type la supernovae. Subsequent 
work on type la supernovae Q, the cosmic microwave 
background and baryon acoustic oscillations @, all 
support the initial observations that the expansion of the 
Universe is accelerating. This late-time acceleration of 
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the Universe is driven by a fluid/field generically called 
dark energy. Very little is known about dark energy. 
Within the context of string theory there is a model for 
the early-time accelerated expansion of the Universe asso- 
ciated with inflation. This string theory motivated model 
for inflation is called Dirac-Born-Infield (DBI) inflation 
and it is driven by the open string sector through 
dynamical Dp-branes. DBI inflation is a special case of 
K- inflation models [12] ■ It was originally thought that 
DBI inflation models would yield large non-Gaussian per- 
turbations which could be used to verify or falsify these 
models and by extension to test string theory [Til 113 ■ 
However, subsequent work has shown that this is may 
not be the case, and that the simplest DBI models are 
effectively indistinguishable from standard field theoretic 
slow roll models inflation [l^ . 

In the present work we examine variants of these DBI 
models as a mechanism, not for the early-time acceler- 
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ation of inflation, but for the observed late-time accel- 
eration. Our DBI scalar field will play the role of dark 
energy. The action for our scalar DBI field is taken to 
have the form found in |l3|. 



to late-time accelerated expansion or not. Of particular 
importance will be stable fixed points which lead to late- 
time accelerated expansion. Finally, we summarize our 
results in section IV. 
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v{<f>) - T{4>) 



(1) 

where we have assumed that the scalar DBI field , 0, is 
spatially homogeneous so that its spatial derivatives can 
be ignored. This is in accord with the fact that dark en- 
ergy seems to be very homogeneously distributed. Note 
that (p has a non-standard kinetic energy term which 
yields a standard kinetic energy term if one expands the 

square root to first order in For a pure AdSs ge- 
ometry with radius R, the warped tension T{(j)) = tc})^ 
is the D3-brane tension with r = l/{gg\) where is 
the string coupling, a' is the inverse string tension and 
A = R^/a' is the 't Hooft coupling in the AdS/CFT 
correspondence. V{4>) is the potential arising from in- 
teractions with Ramond-Ramond fluxes or with other 
sectors Q. Here we take the potential to be quadratic 
V{4)) = m?<lP and the associated D-brane is in the anti- 
de Sitter throat [15]. 

In a spatially flat Fricdmann-Robcrtson- Walker 
(FRW) metric with scale factor, a{t), it can be shown 
that the energy density and the pressure of the 
DBI scalar field are given by 



P0 



^ • ^'^+V{(j}) and 



7+1'' 7 + 1 

where 7 has the form of a Lorentz boost factor 

1 



^ ^ (2) 
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In this paper we analyze this system of a DBI dark en- 
ergy field interacting with dark matter in terms of late- 
time scaling solutions. Such models are different from 
the original work in [l6| which studied the dynamics of 
a DBI field plus a perfect fluid but with no interaction 
between the DBI field and the perfect fiuid. In the course 
of our analysis of this model of DBI dark energy interact- 
ing with dark matter, we find that for certain parameters 
there are late-time attractor solutions or fixed points in 
the phase space of the parameters. 

The organization of the rest of the paper is as follows: 
In section II, we write down two possible phenomeno- 
logical interactions between the DBI dark energy scalar 
field with the dark matter fiuid. In section III, we ana- 
lyze the autonomous equations in terms of the relevant 
variables of DBI dark energy interacting with dark mat- 
ter. We find the fixed points in the phase space flow of 
these variables. We also discuss the stability (i.e. stable 
flxed point, unstable fixed point or saddle fixed point) of 
fixed points and determine if the fixed points correspond 



II. DBI DARK ENERGY SCALAR FIELD 
INTERACTING WITH DARK MATTER 



Cosmological evolution is thought to be largely dom- 
inated by dark energy and dark matter. Dark energy 
gives a gravitationally repulsive effect while dark matter 
is gravitationally attractive. Usually there is no interac- 
tion between these two components such as in the model 
in [isj where graded Lie algebras were used to give a uni- 
fied theory with both dark energy and dark matter, but 
without any interaction between these two components. 

There has been work such as the two measure cos- 
mological model of [l^ where there is some interaction 
between the dark energy and dark matter components 
of the model. However, since the gravitational effects 
of dark energy and dark matter are opposite (i.e. grav- 
itational repulsion versus gravitational attraction) and 
since dark energy appears to be very homogeneously dis- 
tributed, while dark matter clumps around ordinary mat- 
ter, one expects that any interaction between these two 
dark components of the Universe would be weak. In 
this paper we will consider models where there is an in- 
teraction between dark energy and dark matter. The 
dark energy component will come from the DBI scalar 
field in ^ with an energy density of and pressure 
of and the dark matter component will come from 
a fiuid with an equation of state of Wm = Vml Pm = 0. 
Considering a spatially flat Friedman-Robertson- Walker 
background with scale factor, a(i), and allowing for cre- 
ation/annihilation between the DBI scalar field and the 
dark matter fluid at a rate Q we can write down the 
equations for p^ and pm as 



?,H{l + w^)p, 

Pm + 3i?p„ 



-Q 
-Q 



(4) 
(5) 



Here H = a/a is the Hubble rate with derivatives with 
respect to cosmological time, t, being indicated by a dot. 
The DBI scalar field and the dark matter create/decay 
into one another via the common creation/annihilation 
rate iQ. Q represents the interaction between these two 
fields. Although at this point this interaction is generic, 
one can say that if Q > dark energy converts to dark 
matter, _and if Q < dark matter is converted to dark 
energy 

Since there is no fundamental theory which specifies 
a coupling between dark energy and dark matter, our 
coupling models will necessarily be phenomenological, al- 
though one might view some couplings more physical or 
more natural than others. In this paper we consider two 
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types of coupling given by 



Model I Q ^\j^Ppr. 
Model II Q = aHpm 



(6) 
(7) 



where (3 and a are dimensionless constants whose sign 
determines the direction of energy transfer. For positive 
values of the parameters a,/3 > {Q > 0) there is a 
transfer of energy from DBI dark energy to dark matter; 
for negative values of the parameters a, /3 < (Q < 0) 
there is a transfer of energy from dark matter to DBI 
dark energy. 

The interaction given in Model I may be motivated 
within the context of scalar-tensor theories [2l| - |23j where 
similar interaction terms can be found. Generalizations 
of this model allow for /3 = B((j)) and more general forms 
of 1^(0) (e.g. reference [13, Ha])- Interactions of the form 
given by Model II have been considered in [2^ which used 
Q/H — a,nPrn + ^0/50 and in [23| which used Q/H = 

The equation for the rate of change of the Hubble pa- 
rameter is 



H 



1 



[(1 -I- w^) + p„ 



The Hubble parameter is subject to the constraint 



(8) 



(9) 



In this work we use the units such that SttG = 1 where 
G is Newton's gravitational constant. 

We define the fractional density of the DBI dark energy 
and dark matter via 57^ = p^/^iH^ and Q,,ri = Pm/iH^, 
with the condition that £7^ + r2,„ = 1 which comes from 
([9]). The modified Klein-Gordon equation, which follows 
from ([2]) (H]) ([S]) and gives the evolution of the DBI scalar 
field, takes the form 



3g ■ r^ (7 + 2)(7-l) 

72 73 2T (7 + 1)7 



where = AV{(j))/d(l) and = dT{(j))/d(j). Equations 
([5]) ^ (Uni) give a closed system of equations that deter- 
mines the dynamics of the DBI dark energy scalar field, 
0, interacting with dark matter. 

In order to find the fixed points of this systems and 
to study the late-time attractor behavior of these two 
models we introduce the following set of dimensionless 
variables similar to those used in 
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V 



(11) 



With X and y defined in this way we can recover the 
variables for the scalar field model originally proposed in 
[TtI by taking the limit 7 —J- 1. The variable x roughly 
corresponds to the the kinetic energy of the DBI field. 



while y roughly corresponds to the potential energy of 
the DBI field. In addition to these dynamical variables 
X and y we introduce a third variable 7 = 1/7 which is 
connected with the brane tension T{(j)). Taking the in- 
verse of 7 makes the final equations more compact. Thus 
we have exchanged the three variables V{4)) and T(0) 
for x, y and 7. 

In terms of these variables, the Friedmann constraint 
from dl]) can be expressed as 



y 



1 



(12) 



The equation of state of the DBI dark energy is given by 

(13) 



P0 



7X 



y 



y- 



Finally, we introduce two sets of variables related to the 
potential, V{(f), and the brane tension, T((/)). The first 
set, Ai and A2, are defined as 



Ai = - — , 



A. 



T 



(14) 



The second set of variables, Ai and A2, are given by 



Ai = - 



7^-1/2^3/2' 



A, 



yl/2yl/2- 



(15) 



The relationship between the two sets of variables ([141 
and P3|) is given by 



y 7 X 



A2 = 




(16) 



Combining the above definitions, the evolution equations 
for X, y and 7 can be written as the following autonomous 
system: 



Ax 
diV 



diV 

d7 
diV 



Ai2/V3(l-72) 



Q 



2x 



6xH^ 



+ -{7(x2-l)-j/2} 



(17) 



AiyV3(l~7^) , 3;/ , - 2 21 n8^ 
^ + Y|i + 7a;-2;| (18) 



7(1-7') 



3v/(rT^-v/3(l-7)4 f^iy' + ™' 



\2X 



Q 



xV3(l + 7)^^ 



(19) 



where = Ina and 
equations for Ai and A2 



dN 



J_d_ 

H df 



There are also two 



dAz 

dA^ 



-xA2V37(1 + 7) 



TT , 



- 1 



(20) 



(21) 
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Since d\i/dN and dA2/diV can be expressed in terms of 
dx/dN, dy/dN, dj/dN from (HT]) - thus equations 
([20| . (f2T|) are not independent equations and we only 
need to solve three equations for x, y and 7. The time 
evolution equation for H given in ^ can be re-written 
by differentiating the Hubble parameter with respect to 
N which yields 



1 dH 3 r _ n n. 



(22) 



The total effective equation of state for the DBI scalar 
field plus dark matter can be written as 



P4, + P» 

P0 + Pn 



7x^ 



y 



(23) 



In the next section we will investigate the two different 
models, given by ([6]) d?]), near critical points {xc, Vc)- 
Near a critical point the scale factor of the FRW space- 
time takes the form 



(24) 



In order to have accelerated expansion (i.e. a > 0) the 
above equation requires that Weff = ^x^ — y^ < — i Com- 
bining ([22]) with ((231), and recalling that 317 = "ff 3t' 
Hubble parameter evolution equation becomes 



_H_ ^ 3(1 + Weff) 
i/2 2 



(25) 



The energy balance equations (jlj and ([5]) for Model I and 
Model II are independent of H when expressed in terms 
of the variables x{N) and y{N), where iV = Ina. Thus 
the Hubble parameter evolution equation, (|25p . is not 
needed for these particular interacting models, and the 
phase space of both models is a two-dimensional phase 
space involving x and y. 



III. CRITICAL POINTS AND STABILITY 
ANALYSIS 



In this section we find the critical or fixed points of 
the autonomous system pT|) - and perform a sta- 
bility analysis of these fixed points. We are looking for 
the late-time attractor structure of this system of a DBI 
scalar field interacting with dark matter via an energy 
exchange given by Q. The fixed-points for pT)) - 
are found by setting dx/dN = dy/dN = dj/dN — 
and solving the resulting three algebraic equations for 
the critical Xc,yc,^c- Additionally, we will focus on the 
case when 7c = (from equation ^ this implies 7 = 00) 
or 7c = 1 (from equation ([3]) this implies 7 = 1). Thus 7 
is constant and the autonomous system reduces to only 
two dynamical variables - x and y. 

After finding the fixed points we study their stability 
with respect to small perturbations, 6x and 6y, about the 
critical points Xc, yc- Explicitly these take the form 



Xc + Sx, y = yc + Sy 



(26) 



Substituting ([26)1 into (ITTl) and ([T8)) . and keeping terms 
up to first order in Sx and Sy, leads to a system of first- 
order differential equations of the form 



d r Sx 
dN \ Sy 



M 



Sx 
Sy 



(27) 



where is a 2 x 2 matrix that depends on Xc and j/c- To 
study the stability around the fixed points one calculates 
the eigenvalues of M. We denote these eigenvalues as ^1 
and fi2 and for every critical point, {xc,yc), there is an 
associated eigenvalue pair (/ii,/i2)- The stability of the 
critical point, {xc,yc), is then determined by its associ- 
ated eigenvalue pair {fii,fi2) in the following way: (i) If 
/ii < and /i2 < then the critical point is stable; (ii) 
a fii > and /i2 > 0, then the critical point is unstable; 
(iii) if /ii < and ^2 > or > and fi2 < then 
one has a saddle point; (iv) If the determinant of the ma- 
trix is negative and the real parts of ^1 and ^2 are 
negative then one has a limit cycle. For our two models. 
Model I and Model II, all the fixed points fall into cases 
(i) , (ii) or (iii) . None of the fixed point we found are limit 
cycles. 



A. Interacting Model I Q 



Recalling that we are taking 7 to be a non-dynamical 
constant with value or 1 our autonomous system of ()17|) 
- (dH) for Model I © becomes 



dx _ Aiy V3(l - 72) /?^27(1 + 7)[1 - a:;^ - y^ 
dN 



2x 



ix 



+ — {7(x2-l)-z/} 



dy_ 

dN 



7^) 



1(1 



~ 2 
- 7a; 



(28) 



y'} (29) 



The dynamics of this autonomous system is determined 
by the parameters /3 and Ai 

The fixed points for ([28l) (|29|) are obtained by setting 
dx/dN = and dy/dN — and solving the resulting 
algebraic equations for Xc and j/c- There are six fixed 
points and these are presented in the first two columns 
of Table m 

We call the fixed points (al)-(a2) "ultra-relativistic" 
fixed points since for them the Lorentz factor, 7, (jSj 
tends to infinity (which implies that 7 = 0). There are 
also an infinite number of "trivial" fixed points for which 
7c = 2/c = and for which Xc is arbitrary within the range 
constrained by < Xc < \/l ~ {^m is the previously 
defined fractional density of dark matter). These unsta- 
ble, "trivial" critical points are shown along the y = 
axis in figure [1] The four other critical points (bl)-(b4) 
listed in Table |T] are "standard" fixed points since for 
these fixed points the Lorentz factor of ([3]) equals 1 (so 
that 7 = 1) and the DBI field will mimic the behavior of 
a canonical scalar field. 
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TABLE I: The fixed points for Model I with Q = Jl/3p„^(j>. 



r ixed. point 






Vc 


7 


o 








Woff 


(al) 


/ Ai(yA2 + 12-Ai) 
V 6 


^/A 

V 


2 + 12-Ai 
2^3 





1 


[a/a 

L V 


12 


[\/a 

L V 


2 + i2_Ail2 
12 


(a2) 


^/ Ai(V>'? + 12-Aij 




J + 12-Ai 





1 




f + 12-Ai]2 




f + 12->,il2 






2v^ 








12 




12 


(bl) 


-1 







1 


1 




1 




1 


(b2) 







1 


1 


1 




-1 




-1 


(b3) 


1 







1 


1 




1 




1 


(b4) 


2/3 
3 







1 


9 




1 




9 



i. stability of the fixed points in Model I 

For each of the six fixed points listed in Table U we 
found the eigenvalues /ii and /i2 of the matrix M in (|27p . 
The results for each point are listed below along with 
whether the point is stable, unstable or a saddle point 

• Point (al): 

Ml - -\ - ^A2 + 12^ 

M2 = i (^-A2 + Ai^A2 + 12- 12^ 
This point is stable for all values of Ai. 

• Point (a2): 

Ml = -\ (ai - V^A? + 12) 

/i2 = 3 (^-A? + Ai^A2 + 12- 12^ 
This point is stable for all values of Ai. 

• Point (bl): 

Ml = 3, M2 = 3 - 2/3 

This point is a saddle point for /3 > | and is unsta- 
ble for ;3 < |. 

• Point (b2): 

yUl = -3, /X2 = -3 

This point is stable for all values of and Ai. 

• Point (b3): 

/ii = 3, M2 = 2/3 + 3 

This point is unstable for /? > — | and is a saddle 
point for (i > — |. 

• Point (b4): 

A*i = i(4/5'-9),M2 = H4/3' + 9) 

This point is unstable for /? > | or /? < — | and is 

a saddle point for -| < /? < |. 

^. Late-time behavior for model I 

In this subsection we investigate the late-time behavior 
of the scale factor a{t). We are interested if a{t) is accel- 



y 




X 



FIG. 1: Model I: The x - y phase plane for the DBI dark 
energy interacting with dark matter. We have taken the pa- 
rameters as /3 = 1, Al = 1, and 7 = 0. There are two stable, 
critical points - a(l) and a(2). Both points lead to accelerated 
expansion. The dotted line is a; = 0, which is the singularity 
value of the non-allowed region. 

crating, which means that the total effective equation of 
state parameter for this model should satisfy < 
The six critical points from |I] are listed in |ll] with the 
conditions under which they are stable, unstable or sad- 
dle points (this is the second column); the conditions (if 
any) under which the critical point point leads to acceler- 
ated expansion (this is the third column); the conditions 
on /3, Al for the critical point to exist (this is the fourth 
column) . 

Table Hn shows that for Model I, the two critical points 
(al) and (a2) are stable for all value of Ai since /ii, /i2 < 
0. Further these two points lead to accelerated expansion 
(i.e. Wcff < —\) if Al < 2. For these reasons these two 
attractors are of interest in explaining the observed late- 
time accelerated expansion of the Universe. The phase 
space flow for these two points is shown in figure [1] 

The point (b2) is a stable critical point (i.e. ^i,/i2 < 
0) for all values of /? and Ai. Additionally (b2) gives 
accelerated expansion for the Universe (i.e. Woff < — ^) 

for all values of Ai and [3. The x — y phase space behavior 
of the critical point (b2) is shown in figures [5] and |31 
The remaining critical points (bl), (b3) and (b4) are of 
less interest phenomenologically since they do not lead 
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TABLE II: The conditions for stability, acceleration and existence of the fixed point for Model I in terms of the parameters /3 
and Ai. 





Q-f Q Vl! 1 




T^^.YI Qi"PTl PP 
J_jA.io utilise 


Call 


S!i"fiHlp TinHp foT all vfiliipt; nf Ai 


Ai < 2 


all B Ai 


(a2) 


Stable node for all values of Ai 


Ai < 2 


all p, Ai 


(bl) 


Saddle point for ^ > | 


No 


/3 > 




Unstable node for /3 < § 






(b2) 


Stable node for all values of /3, Ai 


Yes 


all /3, Ai 


(b3) 


Saddle point for /3 < -| 


No 


all /3, Ai 




Unstable node for /3 > — | 






(b4) 


Saddle point for — | < /3 < | 


No 


^ < 



Unstable node for 3 < — | or /3 > | 



b2 




FIG. 2: Model I: The x — y phase plane for DBI dark energy 
interacting with dark matter with the parameters = 1, Ai = 
1 and 7 = 1. The point (62) is a stable critical point which 
leads to late time accelerated expansion. 



b2 




FIG. 3: Model I: The x ~ y phase plane for DBI dark energy 
interacting with dark matter with the parameters jS — —1, 
Ai = 1 and 7 = 1. The point (62) is a stable critical point 
which leads to late time accelerated expansion. 



B. Interacting Model II Q = aHp, 



The autonomous system for the variables x and y is 
now 



to accelerated expansion. Thus, there are three fixed- 
points - (al), (a2), and (b2) - which are stable and lead 
to late-time accelerated expansion. However, in regard to 
the existence of these fixed points or their stability the 
specific value of /? (the parameter which characterizes the 
coupling between dark energy and dark matter) plays no 
role. All these fixed points would exist and be stable even 
if /? = i.e. even in the absence of coupling between dark 
energy and dark matter. Thus for Model I the overall 
conclusion is that this coupling does not play a significant 
role in the late time behavior of the system. 



da; 
dN 



y 



AiyV3(l-72) a ^2 ,.2i 



3x 



(30) 



dy AiyV3(l-7^) , 3y , - 2 2 

-r-TT = z 1- — •! H- ^a; - y 

dN 



2 +:^{i + ^.^-/pi) 



The fixed points are again obtained by setting 
dx/dN = and dy/dN = in ([201) and dS]). For Model 
II there are eight fixed points - (dl) to (d4) and (el) to 
(e4) - and these are listed in table IIIII 

The fixed points (dl)- (d4) are "ultra-relativistic" since 
the Lorentz factor, 7, in ([3]) tends to infinity and thus 
7 = 0. The other four fixed points (el)-(e4) are "non- 
relativistic" since the Lorentz factor, 7, equals 1 and thus 
7 = 1. For this four "non-relativistic" fixed points the 
DBI field mimics the behavior of a canonical scalar field. 
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TABLE III: The fixed points for Model II 



Fixed point 




yc 


7 









(dl) 

(d2) 
(d3) 

(d4) 
(el) 
(e2) 

(e3) 

(e4) 



Ai(\/IfTT2-Ai) 



Ai(A/Af + 12-Ai) 
6 

-1 
1 



2^3 



^Af + 12-Ai 









1 
1 
1 

H 

3 

H 

3 



[7Af + 12-Ai]^ 
12 



[VAf + 12-Ai]^ 
12 
1 
1 

1 
1 



[7Af + 12-Ai]^ 
12 



12 
1 
1 

M 

3 

H 

3 



1. Stability of the fixed points in Model II 

The stability analysis of the eight fixed points of Model 
II follows the same procedure as for Model I. For each of 
the eight fixed points listed in Table IIIII we found the 
eigenvalues fii and fj,2 of the matrix Ai in (|?7|) . The 
results for each point are given below along with whether 
the point is stable, unstable or a saddle point 

• Point (dl): 

Ml = a 

This point is cither a saddle point (if a < 0) or 
unstable (if a > 0). 

• Point (d2): 

Ml = = a 

This point is cither a saddle point (if a < 0) or 
unstable (if a > 0) . 

• Point (d3): 

/ii-i(^Ai (^^A2 + 12-Ai^ -12^ 

H2 = fXi (^^Jxl + l2 - Ai^ + a - 3 

This point is either stable (if a < 

i(6 + Af) - ^y^ "'-6a+9 ) unstable (if 

a>i(6 + A?)-iyi^^ ) 

• Point (d4); 

Ml = i (^Ai (^^A2 + 12 - Ai^ - 12^ 

H2 = 5A1 (^^Xj + 12 - Ai^ + a - 3 

This point is either stable (if a < 

i(6 + Xl) - ) or unstable (if 

a>i(6 + A?)-iy^i2^ ) 

• Point (el): 

Ml = 3, M2 = a + 3 



This point is either a saddle point (if a < —3) or 
unstable (if a > —3). 

• Point (c2): 

Ml = 3, M2 = a + 3 

This point is either a saddle point (if a < —3) or 
unstable (if a > —3). 

• Point (e3): 

Ml = -a - 3, H2 = 

This point is a stable (if a > 3 ) , a saddle point (if 
— 3 < a < 3) or unstable (if a < —3). 

• Point (g4): 

Ml = -a - 3, H2 = 

This point is a stable (if a > 3 ) , a saddle point (if 
— 3 < a < 3) or unstable (if a < —3). 



2. Late-time behavior of Model II 

In this subsection we move on to the analysis of the 
late-time attractor structure of Model II as given by the 
autonomous system in ([30|) - (|3T|) . The results for the eight 
fixed points of Model II are summarized in tables Hill and 
IIVI The behavior of the dynamics of the DBI scalar field 
interacting with dark matter via Q = aHpm depends on 
the values of the parameters a and Ai. We found that 
there are non-trivial "scaling solutions" where x, y and 7 
are finite constants depending on the model parameters 
a and Ai. 

Tables IIIII and IIVI show that in the interacting Model 
II, for "ultra- relativistic" fixed points - (dl) and (d2) - 
that /ii is always positive, while fi2 can be either positive 
or negative depending on the value of a. In particular 
(dl) and (d2) are saddles point for a < and a unstable 
for a > 0. For the fixed points (d3) and (d4) fj,i is always 
negative while /X2 can be positive or negative depending 
on the value of a and Ai. In particular (d3) and (d4) are 
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TABLE IV: The conditions for stability, acceleration and existence for the eight critical points of Model II. We list the character 
of the fixed points as a function of the parameters a and Ai. 

Fixed point stability Acceleration Existence 

(dl) Saddle point for a < No for a > and a < 

Unstable node for q > 

(d2) Saddle point for a < No for a > and a < 

Unstable node for q > 



(d3) Stable node for q < i(6 + A?) - f y Ai < 2 for a > and a < 

Unstable node for a > i(6 + A?) - \\f^^^^ 
(d4) Stable point for a < i(6 + Af ) - \ P^I^E+l Ai < 2 for a > and a < 



Unstable node for a > i(6 + A?) - |y 
(el) Saddle point for a < —3 No for a > and a < 

Unstable node for a > —3 

(e2) Stable point for a < —3 No for a > and a < 

Unstable node for a > —3 

(e3) Stable node for q > 3 No for a < 

Saddle point for— 3 < a < 3 
Unstable node for a < —3 

(e4) Stable node for a > 3 No for a < 

Saddle point for —3 < a < 3 
Unstable node for a < —3 




FIG. 4: Model II: The phase plane for the DBI dark energy 
with interaction to dark matter corresponding to a = 1 for 
Ai = 1, 7 = 0. The late-time attractors, (d3) and (d4), are 
stable fixed points which lead to accelerated expansion. 



A?) 



1^ I ct2-6Q+9 



stable, fixed points when a < -^(6 

and are unstable points for a > -i (6 + A^) — | 

These fixed points, (d3) and (d4), are the only critical 
points from Model II which can give rise to an acceler- 
ated expansion for the Universe. Accelerated expansion 
will occur for (d3) and (d4) if Ai < 2. Thus the fixed 
points (d3) and (d4) are good candidates for the late- 
time attractor solution under the requirement that the 
parameters a and Ai are such that these points are sta- 
ble and that they meet the conditions for accelerated 
expansion. The x — y phase space flow of the points (dl) 
- (d4) are shown in figure 2] 

For the "non-relativistic" (i.e. 7 = 1) fixed points 



(el) and (e2) is always positive, whereas /Lt2 can be 
either positive or negative depending on the value of a. 
If a < —3 then /12 < and (el)(e2) are saddle points. 
If a > — 3 then /i2 > and (el)(e2) are unstable fixed 
points. For the points (e3) and (e4) and ^2 can be 
either positive or negative depending on the value of a. 
These points are stable for a > 3, a saddle point for 
—3 < a < 3 and an unstable node for a < —3. Thus 
none of the points (el) - (e4) are stable and from table IIVI 
none of these fixed points lead to accelerated expansion. 
Thus all of these points are not phenomenological viable. 

Only points (d3) and (d4) satisfy the criteria of ac- 
celerated expansion (for values of Ai < 2) and stability 
(for certain values of a and Ai). Here, in contrast to the 
case for Model I, the coupling between dark energy and 
dark matter plays a significant role in the stability or not 
of the fixed points. Although fixed points (d3) and (d4) 
would still exist without the coupling between dark en- 
ergy and dark matter characterized by the parameter a, 
the stability of these points does depend crucially on a 
and therefore on the coupling between dark energy and 
dark matter. From table HVl one can see that if one sets 
Ai = 2 (i.e. the maximum value for which one still gets 
accelerated expansion) then a is restricted as 



2a 



— 6a + 9 



< 10 



Solving this gives the restriction that 
0.088 <a< 4.62 . 



(32) 



(33) 



Only for a in this range (and also one needs Ai < 2) 
does one get both a stable critical point and accelerated 
expansion. 



9 



IV. SUMMARY AND CONCLUSIONS 

In this paper we studied the dynamics of dark energy 
(in the form of a Dirac-Born-Infeld scalar field) interact- 
ing with dark matter (in the form of a fluid) in standard, 
flat FRW cosmology. For the scalar potential and the 
brane tension of the DBI scalar field we took power-law 
functions. We analyzed two different models for the in- 
teraction between dark energy and dark matter. Starting 
with the general interaction equation Q and ([5]) we con- 
sidered Q = ^2/3/3pm0 (i.e. Model I) and Q = aHp^ 
(i.e. Model II). For each model of the dark energy-dark 
matter interaction we were interested in fixed point which 
had accelerated expansion (i.e. We// < ~l/3) and for 
which the fixed points were stable. 

For Model I from tables HI andlUone can see that points 
(al), (a2) and (b2) satisfy these two conditions. However, 
none of these points had any dependence on the coupling 
parameter /? between the dark energy and dark matter. 
Thus for Model I there was no effect of adding the cou- 
pling between dark energy and dark matter or a model 
which only had a DBI scalar field. Thus while not ruled 
out. Model I is not really of interest since the coupling 
to dark matter does not lead to any different result from 
simply having a DBI scalar field. 

For Model II we find two stable fixed points which have 
late time acceleration - points (d3) and (d4) - as can be 
seen in tables IIIII and IIVI In addition from table IIVI one 
can see that whether one has accelerated expansion or 
not depends on Ai which from (dU and ([T5|) depends on 
the scalar field potential and the tension, but does not 
depend on the coupling between the DBI field and the 
dark matter fiuid. Thus addition of a coupling between 
dark energy, in the form of a DBI scalar field and the dark 
matter fluid does not appear to contribute to the exis- 
tence of accelerated expansion for points (d3) and (d4) in 
Model II. A similar conclusion was reached for Model I 
where the stable fixed points with accelerated expansion 
did not depend on the parameter /? which was a measure 
of the coupling between dark energy and dark matter for 
this model. However, for Model II one can see, by looking 
at the "stability" column of table HVl that whether or not 
a given fixed point is stable does depend on a and there- 
fore on the coupling between the DBI field and the dark 
matter fluid. In particular in order for the critical points 
(d3) and (d4) to be lead to accelerated expansion and to 



be stable one needs to restrict the coupling parameter, 
a, between the DBI dark energy and dark matter fluid 
via (p3l) . 

Although we have only examined two speciflc types of 
couplings between dark energy (in the guise of a DBI 
scalar field) and dark matter we will tentatively advance 
some general conclusions about such models: 

(i) The existence of fixed points with accelerated ex- 
pansion does not depend on the coupling between the 
dark energy and the dark matter. One can make a gen- 
eral argument to support this conclusion. Dark matter 
will fall off like [a{t)]^^ while if the dark energy is to 
lead to late time acceleration it will act like a cosmolog- 
ical constant which falls off like [a(t)]° — >■ constant. At 
late times the [a(<)]'' behavior will always dominate the 
[a(t)]~'^ behavior. 

(ii) While the coupling to dark matter does not seem 
to play a role in the existence of fixed points with accel- 
erated expansion it can play a role in their stability. For 
example, in Model II the stability of the fixed points (d3) 
and (d4), which had accelerated expansion, depended 
on the dark energy-dark matter coupling parameter a. 
However for Model I none of the fixed points' stability 
depended on [3, the coupling parameter for this Model 
I. Thus it seems whether coupling between dark energy 
and dark matter is important or not to the stability of 
the flxed points depends crucially on the type of coupling 
one chooses. 
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